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The wahs of bubbles in a first-order phase transition can propagate either as det- 
onations, with a velocity larger than the speed of sound, or deflagrations, which are 
subsonic. We calculate the gravitational radiation that is produced by turbulence 
during a phase transition which develops via deflagration bubbles. We take into ac- 
, count the fact that a deflagration wall is preceded by a shock front which distributes 

i the latent heat throughout space and influences other bubbles. We show that turbu- 

^ \ lence can induce peak values of ^gw as high as ~ 10~^. We discuss the possibility of 

Qh' detecting at LISA gravitational waves produced in the electroweak phase transition 

^ ■ with wall velocities ^ 10~^, which favor electroweak baryogenesis. 

OO 

P^. I. INTRODUCTION 

o 

^ \ One of the possible outcomes of a cosmological phase transition is gravitational radiation. 

^ ■ Although it is difficult to detect, gravitational radiation provides a direct probe of the 
phase transition dynamics, since it propagates freely until the present epoch. Studies of the 
^ , gravitational waves due to phase transitions at the electroweak scale and beyond [l| show 
J> \ that the signal may be within the sensitivity range of LISA and the second generation of 



space-based interferometers. 
^r^ . In a first-order phase transition, bubbles of the stable phase nucleate and expand, convert- 

O \ ing the high-temperature phase into the low-temperature one. Gravitational waves (GWs) 
'nI" I are generated either by the collisions of bubbles a a a 

or by the turbulence that is 

^ ■ produced in the plasma due to the motion of bubble walls [3, a H S S S S 1^1- As bub- 
O ■ bles expand, latent heat is released at the phase boundary. Part of this energy raises the 
temperature of the plasma, and another part is converted to bulk motions of the fluid. For 
phase transitions in the early Universe, the Reynolds number is large enough for the wall 
motion to produce turbulence. In general, the amplitude of GWs coming from turbulence 
, dominates over bubble collisions. 

The bubble wall can propagate via two modes, namely, detonation and deflagration. For 
a detonation, the phase transition front moves faster than the speed of sound Cs, whereas a 
deflagration front is subsonic. GWs produced by detonations have been extensively investi- 
gated Baaaaaaa. In this case, the bubble wall velocity Vw depends only on the 
ratio of the latent heat to the total energy density, a Tl[. Since > Cg, no signal precedes 



the detonation front. Hence, the dynamics of a wall is not influenced by other bubbles 
walls (except in the collision regions). Furthermore, the kinetic energy injected into the 
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fluid is concentrated in a thin region behind the wall. This simplifies the integration of the 
kinetic-energy density profile. Besides, the temperature in the supercooled phase outside the 
bubbles decreases due to the expansion of the Universe. As a consequence, the nucleation 
rate F increases exponentially. This justifies modeling it by F = Fqc^*, which corresponds 
to linearizing the time dependence of the exponent. Hence, (3~^ is the only time scale in 
the problem and determines the duration of the phase transition. At ~ and the mean 
bubble separation, d ~ Vw(3~^. Thus, the final result depends only on the parameters a and 
(3. These features simplify the calculation of GWs. 

The case of defiagrations is more difficult. In contrast to detonations, the defiagration 
front velocity depends on the viscosity of the plasma and on the pressure difference between 
phases, which in turn depends on the amount of supercooling. The wall propagates at a 
subsonic velocity and is preceded by a supersonic shock front which affects other bubbles. 
The shock wave distributes the latent heat injected by the wall throughout space, causing 
a reheating of the Universe and bulk motions of the fluid far away from the wall. For these 
reasons, the case of deflagrations has not received as much attention as the detonation case. 
In Ref. j3] an analytic expression is derived for the signal from bubble collisions. However, 
in that work the fluid velocity is left as a free parameter. Deflagrations have also been 
considered numerically in Ref. [5|. 

The spectrum of GWs for the deflagration case thus depends on several parameters, which 
are difficult to estimate in a specific model. As a consequence, the results for detonations 
have often been used to investigate the GW production in different phase transitions (see e.g. 
IH), disregarding the fact that the phase transition may actually proceed via deflagration 
bubbles. Notice, indeed, that this is likely the general case. For instance, for the electroweak 
phase transition, estimations of the wall velocity give in general subsonic values Vw ~ 10~^ — 
10"^ 0, which are favorable for electroweak baryogenesis [l3j. 



The aim of this work is to address the case of deflagration bubbles with a realistic ap- 
proach. Thus, we take into account the fact that the shock waves coming from a bubble wall 
influence other bubbles. In section [IT] we describe the dynamics of deflagration bubbles. One 
important feature is that, as a consequence of reheating, the nucleation rate turns-off and 
the bubble wall decreases during bubble expansion. Therefore, the approximation F oc e^* 
cannot be used in this case. Furthermore, the smaller velocities at the collision time make 
bubble collisions to be suppressed relative to turbulence as a source of GWs. In section HIT] 
we consider the velocity spectrum of the turbulent fluid that arises from the motion of defla- 
gration walls, and in section HVl we follow the approach of Ref. [l] to derive the gravitational 
radiation produced by turbulence in the deflagration case. In section |V] we estimate the 
dependence of the amplitude and peak frequency on the latent heat, the bubble size, and 
the wall velocity. We show that for a strong enough electroweak phase transition, GWs from 
deflagration bubbles may be detected at LISA. Our conclusions are summarized in section 



II. PHASE TRANSITION DYNAMICS 

A phase transition of the Universe occurs when the free energy of a system depends on 
an order parameter cj) (in general a Higgs fleld), which develops a non-zero value. For a 
flrst-order phase transition, there is a temperature range in which the free energy JF (0, T) 
has two minima separated by a barrier. Generally, the high-temperature minimum is = 
0, corresponding to the symmetric phase, and the low-temperature one is 0m (T) ^ 0, 



corresponding to the broken symmetry phase. At the critical temperature T^, the two 
minima = and = (pm have the same free energy. Below this temperature, 0^ (T) 
becomes the global minimum. At a certain temperature T = Tq, the barrier disappears and 
= becomes a maximum of the free energy. 

The energy density p and entropy density s are related to the free energy density by 
p = Ts + JF and s = —dJ-'/dT. We assume that the phase with = is composed of 
radiation and false vacuum, i.e., JF(0,T) = JF_^ (T) = pa — Pr/S, where pA is the zero- 
temperature energy density of the false vacuum, pA = -^(0 = 0,T = 0) = (0 = 0), with 
V{(f)) the zero-temperature effective potential, and pn is the energy density of radiation, 
Pr = fl'*7r^T^/30, where is the number of relativistic degrees of freedom. Thus, for the 
high-temperature phase we have p+ = Pa + Pr, and for the low-temperature phase we have 
p_ = —TT'_ + JF_. The pressure is given by p = —T. Hence, at T = Tc both phases 
have the same pressure. The latent heat / is defined as the energy density discontinuity 
/ = Ap = p+(Tc) — p_(Tc). Thus, / = TcAs = — TcAjF'. We note that in previous works on 
GW generation the latent heat is sometimes confused with the false vacuum energy density 
(or even with the free energy density). In fact, notice that in general / can be anywhere in 
the range < / < pa + Pr(Tc). On the other hand, it can be shown that Pa is bounded by 
< Pa < pr{Tc)/3 (assuming that the vacuum energy vanishes at T = 0) jl4| . 



The nucleation and growth of bubbles has been extensively investigated (see e.g. [ll|, [12, 
13, 14, 15, 16, [13, 18, 13, 2Gi_21]). According to the standard picture, bubbles of the stable 



phase nucleate with a rate pi 

r ^ T^e-^^/^, (1) 

where (T) is the three-dimensional instanton action, which coincides with the free energy 
of a critical bubble. The nucleation rate is extremely sensitive to the temperature in the 
range Tq < T < Tc. At the critical temperature the radius Rc of the critical bubble becomes 
infinite, so 5*3 = oo and F = 0. In contrast, at T = Tq the radius vanishes, so 5*3 = and 
r ~ T^, which is an extremely large rate in comparison to ~ (T^/Mp)'^. Thus, bubbles 
begin to nucleate at an intermediate temperature Tr for which F ^ H'^. To get an idea of 
the dependence of F on T, consider the thin-wall approximation, which is valid near the 
critical temperature. In this case we can write 5*3 = —AJ-'AttR^/S + aAnRl, where a is the 
surface tension of the bubble wall. Then we obtain the critical radius Rc = 2a/ A and the 
action S3 = 167rcr^/3(AjF)^. For T ^ Tc we can also approximate the free energy difference 
by AjF ^ l(Tc — T)/Tc. Thus, we obtain 

T ^ 3/2(Te-T)2' ^ 

Once a bubble is nucleated, it begins to grow. The bubble radius rapidly becomes much 
larger than Rc. Due to the viscosity of the plasma, the bubble walls immediately reach a 
terminal velocity v^j which is determined by the pressure difference Ap = — AjF and the 
friction with the surrounding particles. It can be approximated by Vyj{T) = AjF(T)/?7, 
where is a friction coefficient. Using again the linear approximation for AJF we obtain 

v^^liTc-T)/r]Tc. (3) 



For strong phase transitions, the latent heat I and the amount of supercooling [Tc — Tp) /T^. 
will be considerable. Notice, however, that the velocity depends also on the viscosity and 
can be subsonic even in this case. 



We shall assume that the wall propagates as a deflagration front. Therefore, the wall 
velocity is lower than the speed of sound in the relativistic plasma, Cg = a/1/3. In this 
case, a shock front precedes the wall with a velocity v^h > Cg. Consequently, the shock 
front of a bubble soon influences other bubbles. For <^ c^, the latent heat is transmitted 
away from the wall and is quickly distributed throughout space. This effect can be taken 
into account by considering a homogeneous reheating of the plasma as latent heat is being 
injected [18|. Even with this simple approximation, the evolution of the phase transition 
must be computed numerically. Therefore, it is not straightforward to find relations between 
thermodynamic parameters such as / or a and the quantities that characterize the dynamics, 
such as the duration of the phase transition or the bubble number density. 

Nevertheless, as we shall see, the general features of the dynamics of slow bubble walls 
provide a good deal of information for the computation of gravitational waves. As a con- 
sequence of reheating, the free energy difference AjF(T) decreases and, according to Eq. 
([3]), the bubble expansion slows down. This effect will be important if the latent heat is 
comparable to the difference 5p = Pr(Tc) — pr{Ty). Thus, if / ^ (5p, the temperature 
will get very close to the critical one^. Then, the velocity will decrease significantly and a 
long phase coexistence stage will take place before the transition completes [3, 19, 20]. In 
any case, due to the exponential dependence of the nucleation rate on temperature, bubble 
nucleation turns off as soon as temperature begins to raise [l^, 2l[. Indeed, notice that the 
temperature T{t) has a minimum, which separates the supercooling and reheating stages. 
According to Eqs. ([T]) and (E]), the nucleation rate has a maximum at that time, and due to 
the extreme dependence of F on Tc — T, the maximum must be a sharp peak. Hence, the 
minimum temperature gives the nucleation temperature Tp. Most bubbles are created in a 
very short time interval 5tr around the time tr corresponding to this temperature. 

For t > t-p, the number density of bubbles rib remains constant. The final size of bubbles 

— 1/3 

is given by ~ . Most bubbles begin to expand at t ~ tr with velocity Vy, = Vi. During 
reheating the wall velocity decreases, and by the time bubbles percolate f „, will be in general 
much smaller than the initial velocity fj. As a consequence, the GW signal from bubble 
collisions will be too low. Therefore, we will consider only GWs from turbulence. The main 
turbulence will be generated during the reheating stage, when is still close to Vi. The 
frequency and amplitude of the GWs will depend on the dynamics of the phase transition, 
which is involved. However, as we shall see, the result will be essentially determined by a 
few parameters. 



III. TURBULENCE FROM DEFLAGRATION BUBBLES 



Turbulence from stirring in cosmological phase transitions has been extensively studied 
js], 0, 0, 0, M, 10|. When turbulence is fully developed, a cascade of energy is established 
from larger to smaller length scales, as eddies of each size break into smaller ones. We define 
the energy dissipation rate per unit enthalpy for a given momentum scale k = 27t/L, 



_ 1 C^Pturb 

w at 



(4) 



This is quite generally the case. For strong phase transitions, we have I ^ p ^ dp. For weaker phase 
transitions, the latent heat is smaller, but so is dp. 



where w = p + p is the enthalpy density, p is the total energy density of the fluid, p is the 
pressure, and pturb = w{v'^)/2 is the kinetic energy density of turbulence. Since the phase 
transition occurs in the radiation dominated epoch, we will make for simplicity the usual 
assumption^ p ~ pr. Hence, p ^ Pr = and w ~ wr = 4/3p/j. If the process is 

stationary, Ek gives the rate at which turbulent energy is received at the scale k from higher 
length scales (and transferred to smaller length scales). If the external source stirs the fluid 
at a single scale ks = 27i/Ls, we have a constant rate Ek = £ for scales L < Ls \7]. In this 
case the turbulent energy in the cascade is characterized by the KolmogorofF spectrum 

E{k)^-^P^^E'/'k-'/'. (5) 
w dk 

The cascade stops at the damping scale ko = 27t/Lo, with <^ Ls, at which the fluid 
viscosity dissipates the injected energy into heat. 

As usual, we assume that the fluid is incompressible and statistically isotropic and homo- 
geneous. Then, the Fourier transform of the velocity has the two point correlation function 

{v,{k)v*{q)) = {27rf5' (k - q) (<5,, - k,k,) P{k), (6) 

where the angular brackets mean a statistical average. Notice that (v2(x)) gives the kinetic 
energy density per unit enthalpy density of the fluid. Therefore, the velocity spectrum P{k) 
can be related to the energy spectrum E{k). The relation is E{k) = /c^P(/c)/27r^ 2^- From 
Eq. ([5]), we have 

P{k) ^ 7r2£2/3A;-ii/3 (7) 

in the inertial range ks <^ k <^ ku- The power spectrum P{k) beyond this range was 
obtained in Ref. [Sj from an ansatz for the real-space correlation function. We have 

( ^L^k"^ for k<^ks 

P{k) ^ 27r{v') I §V3T (I) Lf"k-''l' for ks « k « kn (8) 
[ otherwise, 

and the normalization of the spectrum is given by 

{v') ^ {ELsf , (9) 

which corresponds to the fluid velocity on the largest scale Ls (see below). 

For stationary turbulence, the energy dissipation rate in Eq. (jl]) must equal the power 
that is injected by the source. Thus, for the case of expanding bubbles we have 

E = —Kl — , (10) 

w dt ^ ^ 

where fb is the fraction of volume occupied by bubbles and k is an efficiency factor which 
quantifies the fraction of latent heat that goes into kinetic energy of the fiuid (a fraction 1 — k 



^ In fact, the energy density at the beginning of the phase transition is p+ = pn + pa- For a thermal phase 
transition the false- vacuum energy is bounded by pA < Pr/S and is in general pA ^ Pr [1J|- Part of p+ 
is liberated as latent heat, so that at the end of the transition we have p_ = p+ — I. 



goes into thermal energy and causes reheating of the plasma). The shock wave in front of 
the wall sets the fluid moving outward with a velocity Vf. Hence, the injected kinetic-energy 
density is 

As the phase transition front moves a distance the energy released is proportional to 
IVyjSt. Assuming for simplicity a constant fluid velocity up to a distance ~ c^t, the kinetic 
energy injected in the time 5t will be proportional to PkinCs<5t- Thus, we have 

K = (12) 

We still need to determine the fluid velocity Vf appearing in Eq. (fTT]) . The conservation 
of energy and momentum d^jT'^" = can be applied to the wall discontinuity to obtain 
relations for the quantities on both sides of the wall. Assurning a stationary solution and 



non-relativistic velocities, one obtains the equation (see e.g. [17|,l22|) 



w+v^ = (13) 

where w is the enthalpy, v is the velocity of the fluid in the rest frame of the phase transition 
front, and the plus and minus signs stand for the high- and low-temperature phase regions, 
respectively. In the rest frame of the center of the bubble, the fluid inside the bubble is at 
rest, the wall velocity is v^i = — f_, and the fluid velocity in front of the wall is Vf = v+ —V-. 
Using Eq. (fT3i) we obtain 

w+-w_ Ap + Ap 
Vf = ~ f«„ (14) 

W+ Wr 

where Ap and Ap are the energy-density and pressure differences across the wall. Notice 
that, even if the temperature is homogeneous, p (T) and p{T) are different in each phase. 
Ap{T) and Ap{T) depend on the amount of supercooling. For simplicity, we will assume 
that the temperature remains close to the critical one, so that Ap ^ and the energy density 
discontinuity is given by the latent heat, Ap (Tc) = I. Therefore, we have 

Vf^{l/wR)v^. (15) 

From Eqs. (HH), ([T2]), and dUl) we obtain 

i^=^{^/wr)v^Cs (16) 

Thus, Eq. (fTOl) gives 

e = {Kl/w)h = {l/2)v^csa^h, (17) 

where we have defined the ratio of the latent heat to the enthalpy density, a = 1/wr. 

For the validity of the Kolmogoroff spectrum, it is important that the energy is injected 
at a single scale Ls- For the kind of phase transitions we are interested in, all the bubbles 
are formed in a short interval 6tr around the time tr- Hence, the number density of bubbles 
rif, is set at t ~ tr- For t > tr, the number of bubbles remains constant, and the wall velocity 

(T(t)) is the same for all bubbles. The bubble radius is R{t) = J^^Vwit')dt' , and the 



distribution of sizes has a small width 5R ~ Vw{tT)^tv ^ R- At the moment of collision the 

— 1/3 

bubble size is given by the distance between bubble centers, d ~ 

In fact, turbulence begins as soon as the s/ioc/c fronts collide. The first shocks are emitted 
at t ~ tr and collide when they reach the size after a time ~ d/cg. This time is much 
smaller than the total duration At of the phase transition. Indeed, for deflagrations the 
wall velocity is < Cg. In addition, as explained in section [Tll decreases significantly 
from its initial value Vi. Therefore, we have At ^ d/vi > d/cg. Moreover, as mentioned in 
section [U in the general case the initial velocity may be already -C c^, as suggested by the 
electroweak case. When shock fronts collide, they lose their spherical symmetry. Bubbles 
continue expanding, and their walls continue injecting energy into the fluid until all space 
is filled up. Thus, turbulence begins very soon, at t ~ tr, and remains until the phase 
transition is complete. As bubbles expand, the latent heat is taken away at the speed of 
sound, always stirring the fluid at the same scale Ls ~ d. Hence, it is reasonable to assume 
a Kolmogoroff spectrum with ^5 ~ 2Ti/d. Notice that in the detonation case, in contrast, 
there will be much more small bubbles (with L <^ d) than large ones (with L d), due 
to the exponentially increasing nucleation rate F = Foe'^*. Thus, the small bubbles will in 
principle modify the Kolmogoroff spectrum in the detonation case. 

It is important to estimate the turnover time scale of an eddy, tl ~ L/vl, where vl 
is the characteristic fluid velocity on a length scale L. The rate at which an eddy breaks 
into smaller ones is usually assumed to be roughly ~ r^^. Thus, the rate at which energy is 
transferred in the cascade is e ~ ^l^l^- Therefore, we have vl ~ {eLY^^ and ~ e~^^^L'^^^. 
The turnover frequency is given by Uk = uJs{k/ksf^^ , with ujs = vls/Ls- 

Both Vl and tl increase with L. The maximum length scale is the characteristic scale of 
the source Ls ~ d. In general, d is well inside the horizon 14, 18, 2l[. Hence, the time scale 
for the establishment of a cascade, which is on the order of the maximum turnover time 
Tfi, will be at most on the order of the Hubble time. This justifies neglecting the expansion 
of the Universe in the description of turbulence. On the other hand, it is important to 
compare the characteristic turnover time with the duration of the phase transition. For the 
case of detonations, it has been shown that is always larger than At (see e.g. [g], @|). This 
result is more general. Indeed, since ~ d/vd, At ~ d/v^, and the fluid velocity is always 
smaller than the wall velocity, we have Td > At. Since the stirring source lasts less than r^, 
turbulence is not stationary. Nevertheless, the cascade of energy develops, and it was argued 
[sl that, for the generation of gravitational waves, we can assume stationary turbulence with 
a duration ~ = r^. 



IV. RELIC GRAVITATIONAL WAVES 

We aim to calculate the spectrum we would observe today for GWs originated in a 
cosmological phase transition at time t = and temperature T = T^. As explained in Sec. 
Hlj in the deflagration bubble scenario the wall velocity is in general considerably smaller at 
the time of percolation than at the beginning of bubble expansion. Therefore, we expect 
a weak signal from bubble collisions. In view of that, we shall only consider the source 
provided by turbulence. The calculation of the GW spectrum from primordial turbulence 
has been improved in the last years a a. Here we will use the results of Ref. M (see the 
discussion in Sec. I VI p . In this section we briefly review the derivation and write down the 
results in terms of variables that are suitable for the deflagration bubble scenario. 

For a stochastic background of gravitational waves, the spectrum is characterized by the 



quantity [23| 

now if) = -p^, (18) 

Pc d log f 

where paw is the energy density of the GWs, / is the frequency, and pc is the critical 
energy density today, pc = SHq/ShG. Thus, the gravitational wave energy density per unit 
logarithmic frequency is defined by the relation 

f df dpGw (.^s 

= j '''' 

For a stochastic background, dpcw I d\og f is given by the ensemble average of the Fourier 
amplitudes of the tensor metric perturbation hij. The energy density of gravitational waves 
is 

^ / .^ {dthij{x,t)dthij{:>c, t)) 
IdttG 

where the brackets denote the ensemble average. The source of hij{'x,t) is the transverse 
and traceless piece of the stress-energy tensor, which for a turbulent plasma is given by 
Tij (x, t) = wvi (x, t) Vj (x, t) . Thus, the source for tensor perturbations on each mode k is 
the anisotropic stress Ilij{k), which involves a convolution / d^qvi{q,t)vj(k — c[,t) (see [S] for 
details). Therefore, the energy density spectrum, which involves the average (njj(k)n*j(q)), 
can be related to the velocity spectrum P{k) by means of Eq. ([6]). 

Notice that the velocity correlation function (I6l) does not oscillate in time. However, the 
gravitational radiation is produced by the turbulent eddies, which have a turnover frequency 

= vl/L, with k = 2tt/L. The oscillatory behavior of the source, which is relevant for 
the generation of GWs, is lost in the statistical average. One can account for the turnover 
frequency by replacing fi(k) — > fj(k)e*'^'=* for ks < k < k^. As a consequence, the source for 
hij(k, t) can be modeled as e*^'^*njj(k)9(t — tin)0(ifin — t), where a; = 0;^ for > ks, u = us 
for k < ks, and the Heaviside functions limit the source to the interval tin < t < t&n- We set 
tfin — tin = tls and tfin = t*. For t > t:^, the gravitational wave propagates freely, with the 
dispersion relation uo = k. Hence, the frequency / in Eq. ( |T8l) is given by the wave number 
k of the source. The amplitude of the wave is proportional to Hjj(k) and depends on the 
frequency u of the source. 

The spectrum of GWs is obtained from the quantity {hijCk, t)/ijj(q, t)). It is proportional 
to the dimensionless function 

(fc-2^)^L| + (fc + 2^)^L| 

and to (Hjj(k)H*j (q)). The latter is a four-point spectral function of the velocity, and must 
be reduced in order to use the two-point function ([U]). This is usually done by using Wick's 
theorem, although the velocity field is not Gaussian. We have 

(H,,(k)H^(q)) = ^^^(k - q) 1 d^pP(p)P(|k - p|)(l + + P'), (22) 



where 7 = k-p, j3 = k-k — p. Analytical approximations exist for the integral in Eq. 
(I22I) for large and small scales It can then be evaluated using Eqs. (IHl) and IQ. 



Assuming radiation domination at T = T*, the enthalpy w in Eq. fl22l) can be approximated 
by wr = ipfi/S, and the Hubble rate is given by = SnGpn/S. Finally, at t = t^^ we have^ 

. . f Y3M^s,k){kLsY for k < kg 

PGw{k,Q = (^^{H^Ls)\eLsY'^ I k){kLs)-^'^ for ks < k < ko (23) 



47r 



otherwise. 



The function A{uj,k) gives a different spectral dependence whether the largest eddy ve- 
locity is below or above 1/2, but the GW spectrum always peaks at = ks- For 
deflagrations the velocity of the fluid is smaller than the bubble wall velocity, which is 
smaller than that of sound (c^ ~ 0.58). Thus, only in the limit in which both ~ and 
Vw ~ Cs, we will have > 1/2. We shall deal with wall velocities Vw ^0.1. Therefore, we 
only consider the case < 1/2, for which the energy spectrum can be approximated by 

r {k/ksY/vl^ for k < 2vLsks 
n (h - P^^(H T ^^rl ^4/3 ) ^k/kg) for 2vLsks <k<ks , . 

I otherwise. 

The GWs generated at time t^, redshift due to the expansion of the Universe. The energy 
density scales like a~^, and the frequency like a~^. Therefore, the spectrum today is given 
by Eq. (IMl) . with pnit^) replaced with p_R(to) and the wave numbers replaced with the 
corresponding frequencies. Hence, 

10/3 4/3 \ f < 

^nito) fLs_\ I 4(///,) for 2/^ < / < /, 



- \w) [hi] Kflflr"' for U<f< fJ ^'^^ 

otherwise, 

where flR(to) = Pr/ Pc ~ 4.6 x 10~^ [25^], and fs, fp, fn are the redshifted frequencies 
corresponding respectively to the frequency of the largest eddies fs* = vl^L^^, the peak 
frequency /p* = Lg^, and the dissipation frequency fn* = L'^ . A frequency f^^ redshifted to 
today is given by /o = f^a^^/ao. The ratio of the scale factor at to the scale factor today 
is 

a* „ /100\^/^ lOOGel/ 



- ^ 8 X 10-16 ( — ) 1:2:^, (26) 



It is useful to express /o in terms of f*/H^. The Hubble rate H^, = y/^-nGpR/?) can be 
written as ^ 

H,^2^ WHz ( ( ^ ^ , (27) 

VlOO/ \100GeV J ' ^ ^ 

where we have used the relations G = Mp^, with Mpi = 1.22 x lO^^GeV^, and IGeV ^ 
1.5 X lO^'^Hz. Therefore, the frequency today is given by 

/o = 1.6 X IQ-^'Hz ^ ( ^) A. (28) 

■'^ imGeV VlOO/ ^ ' 



^ In Rcf. [8] the calculations were done using conformal time ry and comoving variables k, x, L, etc. Di- 
mensionless combinations such as HL are readily translated to physical variables, since the scale factor a 
cancels out. 



V. GRAVITATIONAL RADIATION FROM THE PHASE TRANSITION 



As we have seen, the scale Ls is given by the bubble separation. Therefore, the peak 
frequency is 

L = 1.6 X IQ-^mHz ( ^) ^ (29) 

The distance d depends on the dynamics of the phase transition. It can vary from values 
d/H~^ ~ 10~^ for weakly first-order phase transitions, to values d/H^^ ~ 10^^ for strongly 
first-order phase transitions (see e.g. [3, 3, 2^)- For T,, > 100Ge\^, we see that millihertz 



frequencies (corresponding to the peak sensitivity of LISA) are obtained for relatively large 
values of the bubble size, d/H~^ > 10~^ (for ~ 100). 

Setting Ls = d in Eq. fl25l) we obtain, for the maximum of the spectrum today, 

/ J \ 10/3 / \ 4/3 

^ic.l,.. = 3xlO-(^) (^) . (30) 

According to Eq. (fT7|) . the rate e at which energy is injected into the fluid is proportional 
to Vwfb- The speed of bubble expansion may vary considerably due to reheating during 
the phase transition. Nevertheless, we have seen that the turbulence which sources the 
gravitational radiation lasts longer than the phase transition. Hence, the generation of GWs 
is not affected by the details of the time dependence of e, and we can use in Eq. flHUl) a mean 
value e which involves a time average of Vu,{t) fb{t) . A precise calculation of the parameters 
d and e requires a numerical computation of the phase transition. We will address such 



computation elsewhere [26] . Below, we find an approximation for as a function of f ^ and 
d. The bubble separation d also depends on the wall velocity. The larger the initial velocity 
Vi, the quicker the reheating and the sooner the turn-off of the nucleation rate. Roughly, 
we have c? ~ oc Vi [2o| . Unfortunately, the number density of bubbles nb is extremely 

sensitive to the dynamics of the phase transition. This prevents any sensible analytical 
approximation for d as a function of the parameters of the model. 

The bubble expansion rate depends on the wall velocity and the bubble size. Roughly, 
fb ~ Fd'^Vw/d^ = Fv^/d, where F is a geometrical factor (for a spherical bubble, F = An), 
and is the average wall velocity. In general, it will not be a good approximation to take 
a plain time average over the entire duration of the phase transition. The bubble expansion 
will slow down significantly due to reheating. Therefore, the turbulence that is generated 
while the wall velocity Vy^ is close to its initial value Vi will give the main contribution to 
flow- Hence, we can use Vi instead of Vy^ in the approximations above. Inserting Eq. ( |T7|1 
in Eq. (|30ll and using ^ Vi, fb^ Fvi/d, we have 




^Giylpeak - 6 X 10-%av.f^F'/' . (31) 



Again, we see that strong phase transitions favor detection at LISA, since the amplitude is 
maximized for large values of d/H~^ (which also give millihertz frequencies). Weak phase 
transitions give smaller values of d/H~^, and also smaller values of Vi and a. Therefore, we 
will have very low values of flcwifp), together with high values of the frequency, fp > mHz. 
In that case, the generated GWs will not be detected by LISA, but may be detected by 
space-based interferometers of second generation, such as BBO. 



Let us consider an electroweak phase transition at T^, ^ lOOGeV, with ^ 100, and 
set F ~ 10. The bubble separation d depends on the extension of the Standard ModeL 
For the frequency to be in the band that LISA is sensitive to, the phase transition should 
be strongly first-order, so that d/H~^ ~ 10~^. Thus, for a wall velocity Vi ~ 10~^ Eq. 
(pTl) gives flaw ~ 10~^^, provided that the latent heat is large enough (i.e., a ~ 1). A 
large latent heat is consistent with a strongly first-order phase transition. This value of 
^Igw is just at the detection threshold of LISA. Notice, however, that this is an order-of- 
magnitude estimate, and the result is sensitive to several parameters. For instance, the 
enthalpy difference w+ — in Eq. ffT^ will give a larger efficiency factor if it is evaluated 
at the supercooling temperature Tr < Tc. For a phase transition at = ITeV, Eq. fl2U|) 
gives the required frequency fp ~ mHz for a value d/H~^ ~ 10^^. In this case, for a ~ 1 
and Vi ~ 10~^ we have flow ~ 10~^. In Fig. [T] we plot the minimum value of a that is 
needed to achieve a peak value flow if p) > 10"^^ with a peak frequency fp = ImHz for a 
phase ransition that takes place at T = T*. 
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FIG. 1: The values of a and T* that give fp = ImHz and ^Gwifp) = 10 for Vi = 0.1 (solid 
line), Vi = 0.05 (dashed line) and Vi = 0.02 (dotted line). 



VI. CONCLUSIONS 

We have considered the production of gravitational waves in a first-order phase transition, 
due to the turbulence that arises from the motion of deflagration fronts. As we have seen, the 
dynamics of the phase transition is completely different from the case of detonations. The 
main difference is the fact that deflagration walls are preceded by shock fronts, which move 
much faster than them and influence other bubbles. This affects both the bubble expansion 
and the generation of GWs. Firstly, as soon as bubbles begin to nucleate and expand, the 
shock fronts collide and the spherical symmetry is lost. Hence, gravitational radiation can 
be emitted before percolation occurs. Secondly, the quick distribution of latent heat causes 
a global reheating. As a consequence, the nucleation rate turns-off after a very short time, 
and the bubble growth slows down. We have taken these facts into account to calculate the 
spectrum of GWs in a realistic scenario for deflagrations. 



We have shown that the usual assumptions for the turbulent fluid, namely, the validity 
of the Kolmogoroff spectrum and the short duration of turbulence in comparison to the 
Hubble time, apply to the case of deflagration bubbles. We have also derived an analytical 
approximation for the efficiency factor K{a, Vw) for the deflagration case. In previous works 
this factor was either calculated numerically or using rough estimations. This analytical 
approximation makes the treatment of the deflagration case simpler, and will be particularly 
helpful for including the GW calculation in numerical computations of the phase transition. 

We have used the results of Ref. jl] for the GW spectrum from primordial turbulence. 
There, it is assumed that the frequency a; of a gravitational wave is determined by the 
wavenumber of the source mode that generates it, i.e., u = k, rather than being determined 
by the characteristic frequency of that mode, ujl = ^tt/ti. As a consequence, the GW 
spectrum inherits the characteristic wavelength of the source, and the peak is at Up ~ 
Lg^. This is correct for a stochastic and statistically homogeneous source of short duration. 
Recently [oj, it was argued that the turbulent source lasts long enough so that it can be 
treated as stationary. In that case, the resulting GW spectrum would be imprinted with the 
characteristic frequency of the source, and we would have Up ~ us- However, the duration 
of turbulence is on the order of the characteristic turnover time ts ~ cu^^ (since the duration 
of the stirring source is shorter). Therefore, the characteristic frequency and the duration 
of the source of GWs are related by ujsTs = 0(l\ and the typical frequency of the GWs is 
given by the typical wave number of the source |27l |. 

The GW spectrum depends on quantities such as Vy^, fb, and d, which require a numerical 
computation of the phase transition for a more accurate evaluation. We will address such 
computation elsewhere [2(|]. Nevertheless, our estimations show that GWs generated in a 
strongly first-order electroweak phase transition might be detected by LISA. The parameter 
values for which this is possible are roughly constrained by T^, > lOOGeV, d > 10~'^H~^, 
Vw ^ 10~^, a > 0.1. Weakly first-order phase transitions correspond in general to smaller 
values of d/H~^, and a. In that case, the spectrum will have a smaller amplitude and 
the characteristic frequency will be away from the peak sensitivity of LISA. 
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